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Abstract 

In the framework of the theory of open systems based on completely positive quantum dynamical 
semigroups, we give a description of the dynamics of entanglement for a system consisting of two 
uncoupled modes interacting with a thermal environment. Using Peres-Simon necessary and sufficient 
criterion of separability for two-mode Gaussian states, we describe the evolution of entanglement in 
terms of the covariance matrix for a Gaussian input state. We determine the asymptotic Gaussian 
maximally entangled mixed states (GMEMS) and their corresponding asymptotic maximal logarithmic 
negativity, which characterizes the degree of entanglement. Using the symplectic eigenvalues of the 
asymptotic covariance matrix, the expressions of von Neumann entropy and mutual information of 
asymptotic GMEMS are obtained. 

Keywords: quantum entanglement, logarithmic negativity, Gaussian states, entropic measures, open 
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1 Introduction 



Quantum entanglement of Gaussian states constitutes a fundamental resource in continuous variable 
(CV) quantum information processing and communication [I]. For the class of Gaussian states there 
exist necessary and sufficient criteria of entanglement OS] and quantitative entanglement measures 
[31[5].The quantum information processing tasks are difficult to implement, due to the fact that any 
realistic quantum system is not isolated and it always interacts with its environment. Quantum coherence 
and entanglement of quantum systems are inevitably influenced during their interaction with the external 
environment. As a result of the irreversible and uncontrollable phenomenon of quantum decoherence, 
the purity and entanglement of quantum states are in most cases degraded. Therefore in order to 
describe realistically quantum information processes it is necessary to take decoherence and dissipation 
into consideration. Decoherence and dynamics of quantum entanglement in CV open systems have been 
intensively studied in the last years [6HT3]. 

When two systems are immersed in an environment, then, in addition to and at the same time with 
the quantum decoherence phenomenon, the environment can also generate a quantum entanglement of 
the two systems and therefore an additional mechanism to correlate them [9|I141[T5]. In this paper we 
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study, in the framework of the theory of open systems based on completely positive quantum dynamical 
semigroups, the dynamics of the CV entanglement of two modes (two identical harmonic oscillators) 
coupled to a common thermal environment. We are interested in discussing the correlation effect of the 
environment, therefore we assume that the two modes are uncoupled, i.e. they do not interact directly. 
The initial state of the subsystem is taken of Gaussian form and the evolution under the quantum 
dynamical semigroup assures the preservation in time of the Gaussian form of the state. In Sec. 2 
we write the Markovian master equation in the Heisenberg representation for two uncoupled harmonic 
oscillators interacting with a general environment and the evolution equation for the covariance matrix 
of the considered system. By using the Peres-Simon criterion of separability for two-mode Gaussian 
states [2l[l6], we investigate in Sec. 3 the dynamics of entanglement for the considered system. We 
determine the asymptotic Gaussian maximally entangled mixed states (GMEMS) and their corresponding 
asymptotic maximal logarithmic negativity, which characterizes the degree of entanglement. During the 
interaction with the environment, any pure quantum state evolves into a mixed state and for quantum 
information processes it is essential to determine the degree of mixedness of quantum states, which can 
be characterized either by the von Neumann entropy or by the linear entropy [T7]. Using the symplectic 
eigenvalues of the asymptotic covariance matrix, we obtain the expressions of von Neumann entropy 
and mutual information of asymptotic GMEMS, which quantifies the total amount of correlations (both 
quantum and classical) contained in these states. A summary is given in section 4. 



2 Time evolution of covariance matrix for two harmonic oscillators 



We study the dynamics of the subsystem composed of two identical non-interacting oscillators in 
weak interaction with a thermal environment. In the axiomatic formalism based on completely positive 
quantum dynamical semigroups, the irreversible time evolution of an open system is described by the 
following general quantum Markovian master equation for an operator A in the Heisenberg representation 
(t denotes Hermitian conjugation) p^[T9]: 



^[H, Ait)] + ^ ^{2V^^AV, - V}V,A - AV}v^). (1) 



dA{t) 

dt f^^"'-y"" ' 2h 



Here, H denotes the Hamiltonian of the open system and the operators Vj,Vj, defined on the Hilbert 
space of H, represent the interaction of the open system with the environment. 

We are interested in the set of Gaussian states, therefore we introduce such quantum dynamical 
semigroups that preserve this set during time evolution of the system. Consequently H is taken to be 
a polynomial of second degree in the coordinates x,y and momenta Px,Py of the oscillators and yj^Vj 
are taken polynomials of first degree in these canonical observables. Then in the linear space spanned by 
coordinates and momenta there exist only four linearly independent operators Vj=i,2,3,4 [20] : 

Vj = axjPx + ayjPy + bxjX + byjy, (2) 

where axj,cLyj,bxj,byj are complex coefficients. The Hamiltonian H of the two uncoupled identical har- 
monic oscillators of mass m and frequency uj is given by 

H = ^ipl+pl) + ^ix' + y'). (3) 
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The evolution given by a dynamical semigroup implies the positivity of the matrix formed by the scalar 
products of the four vectors a^,., b^^, a^, hy whose entries are the components axj, bxj,ayj, byj, respectively. 
We take this matrix of the following form, where all diffusion coefficients D^x, D^p^,... and dissipation 
constant A are real quantities (we put from now on h= 1): 



V 



Dxx 


-Dxp^ - iX/2 


Dxy 


~DxPy 


Dxp, + iX/2 


^PxPx 


~^ypx 


^PxPv 


Dxy 




^yy 


~^ypy ~ ^ 


— Dxpy 


^PxPy 


-Dypy + iX/2 


^PvPy 



\ 



iX/2 



(4) 



It follows that the principal minors of this matrix are positive or zero. From the Cauchy-Schwarz in- 
equality the following relations hold for the coefficients defined in Eq. 



DxX-Dyy 



^PxPx ^PyPy 



D'xy > 0, 



Kpy ^ 0' 



DxxDp^p^ 



^yy^PyPy 



y 

^Xp^ — ^ ) ^XX^PyPy 

\2 

ypy — 4 ' ^yy^PxPx 



Ky > 0' 



D 



ypx 



> 0. 



(5) 



We introduce the following 4x4 bimodal covariance matrix: 



/ (^xxit) CTxp^ii) 0-xy{t) <^xpy{t)\ 



ait) 



XP'. 



a 



xy 
xpy 



it) 



'Xp: 
(^PxPxit) 

(^ypxit) 

'PxPy ^ 



yy 
ypx 
yy\ 



a,iy{t) 



'XPy\ 

'PxPy ^ 

'^ypy^ 



0"„,-n.,(t) 



'ypy 



it) 



cr. 



(6) 



PyPy 



it)j 



The problem of solving the master equation for the operators in Heisenberg representation can be trans- 
formed into a problem of solving first-order in time, coupled linear differential equations for the covariance 
matrix elements. Namely, from Eq. ([T]) we obtain the following system of equations for the quantum 
correlations of the canonical observables, written in matrix form [20j: 

d(T{t) 



dt 



Ya{t) + a{t)Y' +2D, 



(7) 



where 





/ 


-A 


1/m 





\ 




( Dxx 


Dxpx 


Dxy 


Dxpy ^ 


Y = 






-A 










-^XPx 


Dpxpx 


^ypx 


DpxPy 












-A 


1/m 




Dxy 


^ypx 


Dyy 


^ypy 




I 








— mw^ 


-V 




\Dxpy 


^PxPy 


^ypy 


^PyPy / 



(8) 



The time-dependent solution of Eq. ([7]) is given by [2D] 

a{t) = M{t)[a{0) - a{oo)]M^{t) + cj(oo). 



(9) 



where the matrix M{t) = exp(yt) has to fulfill the condition lim^^oo M(t) = 0. In order that this limit 
exists, Y must only have eigenvalues with negative real parts. The values at infinity are obtained from 
the equation 

Ya{oo) + a{oo)Y^ = -2D. (10) 
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3 Dynamics of continuous variable entanglement 



3.1 Time evolution of entanglement 

The characterization of the separabihty of CV states using second-order moments of quadrature operators 
was given in Refs. [2l[3]. A two-mode Gaussian state is separable if and only if the partial transpose of 
its density matrix is non- negative [necessary and sufficient positive partial transpose (PPT) criterion]. 
A two-mode Gaussian state is entirely specified by its covariance matrix ([6]), which is a real, symmetric 
and positive matrix with the block structure 



^(0 = . ' (11) 




where A, B and C are 2x2 Hermitian matrices. A and B denote the symmetric covariance matrices 
for the individual one-mode states, while the matrix C contains the cross-correlations between modes. 
Simon [2] derived the following PPT criterion for bipartite Gaussian CV states: the necessary and 
sufficient criterion of separability is S{t) > 0, where 

= detAdetS+ Q-ldetCl") - T^[AJ C J B J J] - h^Aei A + dei B) (12) 



and J is the 2x2 symplectic matrix 



Since the two oscillators are identical, it is natural to consider environments for which D^x = Dyy^ D^p^ = 
Dypy, Dp^p^ = Dpypy, Dxpy = Dyp^. Then both unimodal covariance matrices are equal, A = B, and the 
entanglement matrix C is symmetric. 

In order to describe the dynamics of entanglement, we used the PPT criterion [2l|T6] according to 
which a state is entangled if and only if the operation of partial transposition does not preserve its 
positivity. We analyzed the time evolution of the Simon function S{t) (jl2p for a thermal environment 
characterized by the temperature T and the diffusion coefficients 

muDxx = ^^^^ = ^ coth Dxp, = 0, m^w^D^^ = Dp^p^, (14) 

corresponding to the case when the asymptotic state is a Gibbs state [19]. We described the generation and 
evolution of entanglement in terms of the covariance matrix for a Gaussian input state and showed that for 
some values of the temperature of environment, the state keeps for all times its initial type: separable or 
entangled. Depending on the values of temperature, mixed diffusion coefficients and dissipation constant, 
entanglement generation, entanglement sudden death or a repeated collapse and revival of entanglement 
can take place [21] . 
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3.2 Logarithmic negativity 

In order to quantify the degree of entanglement of the infinite-dimensional system states of the two oscil- 
lators it is suitable to use the logarithmic negativity. For Gaussian states, the measures of entanglement 
of bipartite systems are based on some invariants constructed from the elements of the covariance ma- 
trix [6l[TU]. Heisenberg uncertainty principle can be expressed as a constraint on the global symplectic 
invariants A = det A + det B + 2 det C (seralian) and det a j2j : 

A<^ + 4deta, (15) 

and the symplectic eigenvalues z/zp, which form the symplectic spectrum of covariance matrix a are 
determined by 



2^1 = A:^ a/A2 -4deta. (16) 

In terms of z^zp relation ()15p takes the form > 1/2. A two-mode Gaussian state is pure if and only if 
deto" = 1/16. Equivalently, the necessary and sufficient criterion for a state to be pure is = = 1/2. 
We notice also the useful relations 

det a = A = v'i + i^l. (17) 

For a Gaussian density operator, the negativity is completely defined by the symplectic spectrum of 
the partial transpose of the covariance matrix and it is given by 

£;7V = max{0,-log2 2P_}, (18) 

where i/- is the smallest of the two symplectic eigenvalues of the partial transpose a of the two-mode 
covariance matrix a : 



2P^ = A =F \/A2 -4deto-. (19) 

Here A is the symplectic invariant, given by A = A — 4 det C = det A + det B — 2 det C. 

Logarithmic negativity determines the strength of entanglement for E]\f{t) > 0. A state is separable 
if and only if 

P„ > 1/2 (20) 

and logarithmic negativity quantifies the violation of inequality ()20p . As expected, the time evolution of 
the logarithmic negativity has a behaviour similar to that one of the Simon function in what concerns 
the characteristics of the state of being separable or entangled 



3.3 Asymptotic entanglement 

Prom Eqs. (llOh and (I14p we obtain the following elements of the asymptotic matrices yl(oo) = B{oo) : 

muja^^{oo) = = Icoth-^, cT^p,(oo) = (21) 

mu! 2 2ki 
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and of the entanglement matrix C(oo) 



(^xpyioo) = (Typ^ioG) = ^^2 ' (23) 



^P.Py (oo) = ^ — ■ (24) 



The mixedness of a quantum state p is characterized by its purity fi = Trp^. For a two- mode Gaussian 
state the purity is given by ^ = l/4\/deto-. The marginal purities pi {i = 1,2) of the reduced states 
in mode i are given by /ii = l/2VdetA and p2 = l/2^/detB. The global and marginal purities range 
from to 1, and they fulfill the constraint fi > piP2, as a direct consequence of Heisenberg uncertainty 
relations. In Ref. |26j the following upper and lower bounds on the invariant A have been obtained in 
terms of global and marginal purities: 



- + ^^^^^<A<min|i^^i±i#-J-,ifl + i,)L (25) 



The invariant A has a direct physical interpretation [26]: at given global and marginal purities, A 
determines the amount of entanglement of the state and the smallest symplectic eigenvalue D- of the 
partially transposed state is strictly monotone in A. Consequently, the entanglement of a Gaussian state 
with fixed global purity p and marginal purities pi, P2 is strictly increasing with decreasing A. According 
to double inequality ()25p . giving lower and upper bounds on A, there exist both maximally and minimally 
entangled Gaussian states. 

We analyze now the existence of the entanglement in the asymptotic regime in the considered sym- 
metric situation A = B and in the particular case D^y = 0. In the limit of long times, we obtain: 

detA = detB = ^, detC7 = -^, deta=(^-^], (26) 
4 A2 V 4 A2 / 



where we used the notations: 



C^T^coth^, d = D,p^, A^ = Lo^ + X\ (27) 



Then we get 



1 = C'2 4^, — = — = Ct, A = 2\/d^=— . (28) 
p pi p2 2/x 



These values saturate the lower bound in inequalities (I25p and this situation entails a maximal entangle- 
ment. Consequently, the corresponding states are Gaussian maximally entangled mixed states (GMEMS). 
They are thermal squeezed states with the squeezing parameter r given by tanh2r = d/ACx- In the pure 
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case (det a = 1/16), these states are equivalent to two- mode squeezed vacua with the squeezing parameter 
determined only by the temperature of the environment : tanh2r = C|, — 1/ 2Ct ■ 
According to Ref. [26], these states are separable in the range 

/X < ^ . (29) 

IJ'l+ fJ'2- 

Then, for a given temperature T, we obtain that the asymptotic final state is separable for the following 
range of positive values of the mixed diffusion coefficient d: 

^<Ct-1. (30) 

We remind that, according to inequalities ([5]), the coefficients have to fulfill also the constraint ACr/2 > d. 
For a given temperature of the environment and for this range of mixed diffusion coefficients, no entan- 
glement can occur for these states. Outside this region, i.e. for a temperature and diffusion coefficient 
satisfying 

Ct-1<^ <Ct + 1, (31) 

they are GMEMS. 

The asymptotic logarithmic negativity has the form 

EM{oo) = -log2(cT-^y (32) 

Outside the separable region, this is the maximum possible value of the logarithmic negativity, attained 
by GMEMS. It depends only on the mixed diffusion coefficient, dissipation constant and temperature, 
and does not depend on the parameters of the initial Gaussian state. 

3.4 Entropy and mutual information 

The degree of mixedness of a quantum state p can be characterized either by von Neumann entropy 
Sv{p) or by the linear entropy Sl{p). For CV systems these quantities are defined as 

Sv{p) = -Tiplnp, Sl{p) = 1 - Trp^ = 1 - p. (33) 

In Ref. [17] it was proven that von Neumann entropy 5'v(o") of an arbitrary two- mode Gaussian state 
with covariance matrix a has the following expression: 

5y(a) = /(/._) + /(/.+), (34) 

with given by Eqs. ([T6]) and 

fix) ^{x + \) ln(x + 1) - (x - ^) ln(x - ^). (35) 

In the symmetric case considered in the previous Subsection, corresponding to a maximal entangle- 
ment, we obtain the following expression of the von Neumann entropy of GMEMS: 

(a) = (2i/M + 1) Hi^Ai + 1)- {2vM - 1) H^M -\), (36) 
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where 



UM = u- = u+ = ^^-j^. (37) 

We see that von Neumann entropy of a two-mode Gaussian state depends on the two invariants A and 
deter, whereas the purity depends only on deto". This imphes that the hierarchy of mixedness estabhshed 
by von Neumann entropy on the set of two-mode Gaussian states differs from that induced by the hnear 
entropy. Namely, there may exist states with a given linear entropy, i.e. with a given det a, but with 
different von Neumann entropies, i.e. with different A. Therefore von Neumann entropy provides a richer 
characterization of the lack of information of a state [17] . 

Knowing the von Neumann entropy one can obtain the mutual information /(p) = Sv{pi) + Sv{p2) — 
Sv{p) (here pi is the reduced density matrix of subsystem i = 1,2), which quantifies the total amount 
of correlations (both quantum and classical) contained in a state J7]. The mutual information /(c) of a 
Gaussian state with covariance matrix a is defined as 



I{a) = Sv{ai) + Sv{a2)-Sv{(7), (38) 

where denotes the reduced one- mode state obtained by tracing over subsystem j ^ i- For an arbitrary 
two-mode Gaussian state the mutual information 1(a) is given by [17j : 

I(ct) = /(a) + f{b) - f{u^) - f{u+), (39) 

where a = \J det ^4, h = Vdet B, and f{x) is given by Eq. (j35]) . 

For a symmetric state with covariance matrix as, i.e. a state with a = b, and for the case corresponding 
to a maximal entanglement, we obtain for the mutual information of GMEMS the following expression: 

I^ias) = (2a + 1) ln(a + ^) - (2a - 1) ln(a - ^) - (2z/m + 1) Hi^m + ^) + (2z/m - 1) ln{i^M - ^), (40) 
where a = Ct/2. 



4 Summary 



In the framework of the theory of open quantum systems based on completely positive quantum 
dynamical semigroups, we investigated the Markovian dynamics of the quantum entanglement for a 
subsystem composed of two noninteracting modes embedded in a thermal environment. By using the 
Peres-Simon necessary and sufficient criterion of separability for two-mode Gaussian states, we have 
described the entanglement in terms of the covariance matrix for Gaussian input states. For a given 
temperature, we calculated the range of mixed diffusion coefficients which determine the existence of 
asymptotic Gaussian maximally entangled mixed states (GMEMS). We obtained also the expression 
of the maximal logarithmic negativity, which characterizes the degree of entanglement of these states. 
Using the symplectic eigenvalues of the asymptotic covariance matrix, we have characterized mixedness 
and total amount of correlations (both quantum and classical) contained in asymptotic GMEMS by 
determining von Neumann entropy and mutual information for these states. 
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